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ABSTRACT 


DNA sequence data obtained from the 96 year old holotype of Marmosa elegans coquim- 
bensis Tate, 1931, support the hypothesis that this nominal taxon is a synonym or subspecies 
of Thylamys elegans (Waterhouse, 1839) and is not conspecific with T: pallidior (Thomas, 1902). 


INTRODUCTION 


The nominal taxon Marmosa elegans coquimbensis was named by Tate (1931) based on a single 
specimen in the Field Museum of Natural History (FMNH) collected on 20 June 1923 by Colin C. 
Sanborn at or near Paiguano (30°02' S, 70°27’ W, ca. 1600 m elev.) in the Coquimbo region of Chile. 
Tates trinomial usage was maintained by most subsequent researchers (e.g., Cabrera, 1958) until 
elegans Waterhouse, 1839, and other species of fat-tailed mouse opossums were transferred to the 
genus Thylamys Gray, 1843, by Gardner and Creighton (1989). Most subsequent authors have 
treated coquimbensis as a subspecies (Palma, 1997) or a junior synonym (Creighton and Gardner, 
2008) of T: elegans, but we, together with S. Jansa (Giarla et al., 2010), transferred coquimbensis from 
the synonymy of T: elegans to that of its sister species T: pallidior (Thomas, 1902) based on morpho- 
logical traits of the holotype (FMNH 22302). If valid, our reidentification of FMNH 22302 as T. 
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pallidior would have been the first record of that species from the Coquimbo region. Recently, 
however, Boric-Bargetto et al. (2016) analyzed cytochrome b sequence data from several specimens 
of Thylamys collected near Paiguano and recovered them as members of a northern haplogroup of 
T. elegans. Although the authors admitted the possibility that T. elegans and T: pallidior might be 
sympatric near Paiguano, they concluded that the haplogroup in question should be called T. e. 
coquimbensis, In order to test their conclusions, we extracted and sequenced DNA from the 96 year 
old holotype of coquimbensis for phylogenetic analysis. 


Materials and Methods 


A single claw was clipped from the dried study skin of FMNH 22302, To avoid contamination 
from exogenous DNA, all pre-PCR laboratory procedures were performed in a biological safety cabi- 
net with UV sterilization in a lab where mammalian DNA is never amplified and in which contami- 
nating mammalian PCR products are unlikely to be present, Moreover, no other didelphid DNA 
samples had ever been present in the building where the work was conducted prior to the experi- 
ment. The sample was soaked in ethanol overnight and then rinsed with water three times. The same 
wash procedure was performed the next day. On the third day, the sample was allowed to air-dry 
before DNA extraction was performed with a modified DNeasy Blood and Tissue Kit (Qiagen) 
protocol. The sample was added to a mixture of 40 ll proteinase K, 160 ul Buffer ATL, and 30 wl 1 
M dithiothreitol and incubated at 56° C in a shaking mixer for 24 hours. The digested sample was 
then lysed by adding 200 ul Buffer AL, mixing, and incubating at 70° C for 10 minutes. DNA was 
precipitated using 200 ul cold ethanol and incubated at 4° C for one hour, The sample was added to 
a QIAquick PCR Purification Kit (Qiagen) spin column, This column was chosen over the standard 
DNeasy spin columns because it is designed to preferentially bind fragments of DNA less than 10 kb. 
The subsequent wash and elution steps followed the usual DNeasy tissue extraction protocol. 

Nested pairs of primers were designed to span 174-359 bp pieces of the mitochondrial cyto- 
chrome b (CYTB) gene in the degraded sample. To ensure that primers would amplify both a rela- 
tive of Thylamys pallidior and T: elegans, primers were designed using a consensus sequence of two 
T. pallidior CYTB sequences (GenBank accession numbers HM583398 and HM583387) and one 
T. elegans whole mitochondrial genome sequence (NC005825). Primers were designed in Geneious 
R9 (Biomatters, Inc.) using the Primer3 algorithm (Untergasser et al. 2012), and only regions that 
exhibited sequence conservation between the two reference species were used (table 1). Seven sets 
of PCR reactions were conducted using seven primer pairs. Each reaction contained 13 ul of GoTaq 
Green Master Mix (Promega), 9 ul of water, 1 jl of each 10 4 M primer solution, and 1 wl of undi- 
luted sample DNA. A negative control reaction, where no DNA was added, was conducted for each 
separate PCR mixture. The reaction mixture was PCR amplified using an initial 2 min. melting 
phase at 95° C; then 40 cycles of 30 s at 95°, 30 s at 55°, and 30 s at 72°; then 5 min at 72°. PCR 
product size was verified on a 1% agarose gel. PCR products were purified using ExoSAP-It (Thermo 
Fisher) and sent to GENEWIZ (South Plainfield, NJ) for Sanger sequencing in both directions. 

For each PCR amplicon, chromatograms were examined, trimmed, and assembled in 
Geneious. To determine whether individual amplicons might be derived from contaminant DNA, 
each was subjected to a standard nucleotide BLAST search (Altschul et al., 1990) against Gen- 
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TABLE 1. Pairs of primers used to sequence overlapping pieces of cytochrome b and associated information 
about PCR amplicons. 


Amplicon Forward Primer Reverse Primer Amplicon GenBank 
number Length Accession of 
Best BLAST 
Hit 
1 Coqui-CY TB-32F; Coqui-CYTB-331R: 300 KP994523.1 
TGGCATGAAAAACCATTGTTGT TCCATTGGCGTGAATGTTTCG 
2 Coqui-CY TB-222F: Coqui-CYTB-511R: 290 KP994527.1 
TCCTAGCCATACATTACACATCAGA = AATTACTGTAGCTCCTCAGAAGG 
3 Coqui-CY TB-445F; Coqui-CY TB-739R: 295 KP994526,1 
AGCTACTGCATTTGTAGGTTATGT GGGGTGGAAAGGAATTTTGTCTG 
4 Coqui-CY TB-457F: Coqui-CYTB-815R: 359 KP994526.1 
TCCTTCTGAGGGGCTACAGT GGGGTGAAGTTATCGGGGTC 
5 Coqui-CYTB-717F; Coqui-CYTB-1016R: 300 KP994522.1 
CAGACAAAATTCCTTTCCACCCC AGCTTCGTTGGTTTGATGTATGA 
6 Coqui-CYTB-901F; Coqui-CYTB-1191R: 291 KP99452651 
TGCCTACGCAATTCTACGATCT TCAAGTATGCCTGCTAAGGGT 
7 Coqui-CYTB-1079F: Coqui-CYTB-1293R: 215 AJ508401.2 


ACCTGAATTGGAGGACAACC 


TGTCCAAGGAAGGAGTTTTCC 


Bank's nonredundant nucleotide database. After ruling out contamination, FMNH 22302 CYTB 
amplicons were assembled to a reference Thylamys pallidior sequence and trimmed to include 
only the CYTB coding region. A representative set of 72 Thylamys CYTB sequences and one 
sequence from the outgroup Lestodelphys halli were downloaded from GenBank (appendix 1). 
This set included sequences from all the Thylamys species recognized by Giarla et al. (2010), 
broad geographic sampling within the Elegans Group, and all the sequences that Boric-Bargetto 
et al. (2016) identified as T: elegans coquimbensis. These sequence data were aligned in Geneious 
using MUSCLE (Edgar 2004). The alignment was partitioned by codon position, and different 
partitioning schemes and nucleotide substitution models were tested in PartitionFinder2 (Lanfear 
et al. 2016). A Bayesian phylogeny was inferred in MrBayes 3.2.6 (Ronquist et al., 2012) using the 
best-fitting partitioning scheme and substitution models. Two simultaneous Markov chain Monte 
Carlo (MCMC) runs were initiated, each for 10 million generations and with four chains. Con- 
vergence of the MCMC runs was assessed in Tracer 1.7.1 (Rambaut et al., 2018) to ensure that 
each parameter estimate had an effective sample size greater than 200. 


RESULTS 


A complete CYTB sequence from FMNH 22302 (GenBank accession MK907779) was assem- 
bled from seven individually amplified and overlapping PCR amplicons derived from high-quality 
Sanger reads. When subjected to individual BLAST searches, each amplicon matched a Thylamys 
elegans sequence in the nonredundant nucleotide database (table 1), ruling out contamination from 
nondidelphids. Once assembled, the CYTB amplicons overlapped each other by 84.5 bp on average 
(range: 20-218 bp), with no conflicting bases. After alignment, PartitionFinder2 identified a best- 
fitting partitioning scheme and associated nucleotide substation models in which each codon posi- 
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Chile, Coquimbo: NK96571 (KP994526) 
°8ILL Chile, Coquimbo:UCK23 (KP994525) 
oa, [¢ Chile, Coquimbo:UCK20 (KP994524) 
Chile, Coguimbo: FMNH22302 (MK907779) 
Chile, Coquimbo: II114 (KP994522) 
Chile, Coquimbo: UCK809 (KP994523) 
Chile, Atacama: UCK21 (KP994521) 
Chile, Coquimbo: NK27583 (HM583376) 
Chile, Coquimbo: NK95622 (KP994527) 
Chile, Atacama: CZZA-UTA252 (KP994520) 
a Chile, Coquimbo: NK96879 (KF164531) 
Chile, Taparaca: SSUCMa519 (KF 164532) 
Chile, Taparaca: SsUCMa520 (KF 164533) 
Chile, Valparaiso: NK96763 (HM583379) 
, Chile, Metropolitana: NK95111 (HM583378) 
Chile, Metropolitana: NK95691 (KF 164528) 
Chile, Metropolitana: NK95436 (KF164526) 
0299 Chile, Valparaiso: UCK807 (KP994529) 
Chile, Valparaiso: NK105407 (KP994528) 
Chile, Santiago: NK95677 (KF164529) 
Chile, Valparaiso: UCK22 (KP994531) 
Chile, Valparaiso: NK27606 (HM583377) 
Chile, Santiago: NK95354 (KF164527) 
Chile, Metropolitana: NK96791 (HM583380) 
Chile, Santiago: NK95971 (KF164530) 
Chile, Maule: NK160945 (KF 164538) 
0.9/4 Chile, Maule: NK160972 (KF 164539) 
I Chile, Maule: NK106178 (KF164537) 
Chile, O'Higgins: NK105928 (KF 164540) 


0.94. 


0.8 


0.97, 


a | Chile, Maule: NK160518 (KF164535) 
0.05 Z Chile, Maule: NK160526 (KF164534) — Thylamys 
substitutions per site 099° L Chile, Maule: NK160466 (KF 164536) elegans 
————___| Thylamys pallidior 
0.99 
<—_| Thylamys sp. Giarla et al. 2010 
0.57 
<—_| Thylamys tatei 


—————]] Thylamys pusillus 
<j Thylamys macrurus 
<——] Thylamys sponsorius 


=<———] Thylamys venustus 


<] Thylamys velutinus 
Thylamys karimii 


FIG. 1. Bayesian phylogenetic tree of Thylamys cytochrome b sequences. Numbers at nodes indicate posterior 
probabilities (PP). Filled circles at nodes denote PPs equal to 1.0. Unmarked nodes received PPs less than 0.5. 
Within T. elegans, tips are labeled with country, region, specimen identifier, and, in parentheses, a GenBank 
accession number. The holotype of Marmosa elegans coquimbensis ‘Tate, 1931, is in boldface type. For other 
species, tips of the phylogeny are collapsed and the outgroup is not shown. See appendix 1 for a full list of 
sequences included in the phylogeny. A full tree file corresponding to this topology is available on TreeBase 


(doi: http://purl.org/phylo/treebase/phylows/study/TB2:S25505). 
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tion was allowed its own substitution model (position 1: GTR+G; position 2: GTR+]; position 3: 
GTR+I+G). The resulting phylogeny unambiguously recovered the holotype of coquimbensis as 
nested within T elegans as that species is currently recognized (fig. 1). 


DISCUSSION 


Our results unequivocally support Boric-Bargetto et al’s (2016) conclusion that coquim- 
bensis is more closely related to Thylamys elegans than to T. pallidior, although the gene tree 
we recovered does not address the issue as to whether coquimbensis should be treated as a 
subspecies or as a synonym of the former species. Indeed, the trinomial classification of what 
is currently considered T. elegans is complicated by several divergent but unnamed mtDNA 
lineages, one consisting of two sequences from Tarapaca (SSUCMa519, SSUCMa520; the “Loa 
lineage” of Boric-Bargetto et al., 2016), another of a single sequence from Valparaiso (NK96763), 
and a third of seven sequences from Maule and O’Higgins (NK105928, NK106178, NK160466, 
NK160518, NK160526, NK160945, NK160972), Of these, only the “Loa” haplogroup from 
Tarapaca was analyzed by Boric-Bargetto et al. (2016), although the Maule/O’Higgins and 
Tarapaca sequences had previously been analyzed by Palma et al. (2014). 

Our results also raise questions about phenotypic character variation within Thylamys elegans. 
The type of coquimbensis (FMNH 22302) is—as previously reported by us (Giarla et al., 2010: 45)— 
phenotypically very similar to T. pallidior, notably in pelage coloration (much paler dorsally than 
typical T: elegans, and with almost entirely self-white ventral fur), hind-foot length, and bullar infla- 
tion. Apparently, the Paiguano topotypes examined by Boric-Bargetto et al. (2016) are also pale-furred, 
but other specimens that they refer to T. e, coquimbensis (e.g., NK 27583, one of a series of specimens 
that we examined at the Museum of Southwestern Biology in Albuquerque, NM) are much darker, 
so this trait is evidently not diagnostic of the haplogroup that they associate with this trinomen. 

The other traits of FMNH 22302 that led us to synonymize coquimbensis with Thylamys pal- 
lidior (smaller feet, more inflated bullae) were not mentioned by Boric-Bargetto et al. (2016), wha 
emphasized the diagnostic value of nasal morphology, in which their Paiguano specimens were 
said to resemble T: elegans rather than T. pallidior. We had previously evaluated this character 
(which was also mentioned by Solari, 2003) as a diagnostic criterion for our revision (Giarla et 
al., 2010), but found too much intraspecific variation for it to be useful, and we invite readers to 
inspect Boric-Bargetto et al’s (2016: fig. 2) illustration to decide this issue for themselves. 

In summary, the mtDNA sequence data at hand clearly support Boric-Bargetto et al’s 
(2016) decision that the nominal taxon coquimbensis should be associated with Thylamys ele- 
gans rather than with T! pallidior, but there remain unresolved issues of trinomial nomenclature 
and geographic variation that should be addressed in any future revision of these taxa. 
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APPENDIX 1 


CYTOCHROME B SEQUENCES DOWNLOADED FROM GENBANK AND USED IN THE 
PHYLOGENETIC ANALYSIS 


Species Museum Code Specimen Identifier Genbank Accession Number 
Lestodelphys halli PNG1399 KF164579 
Thylamys elegans CZZA-UTA252 KP994520 
Thylamys elegans Ui14 KP994522 
Thylamys elegans MSB207075 NK105407 KP994528 
Thylamys elegans MSB153832 NK105928 KF164540 
Thylamys elegans MSB 207538 NK106178 KF164537 
Thylamys elegans NK160466 KF164536 
Thylamys elegans MSB 230411 NK160518 KF164535 
Thylamys elegans MSB 230254 NK160526 KF164534 
Thylamys elegans NK160945 KF164538 
Thylamys elegans NK160972 KF164539 
Thylamys elegans MSB 87097 NK27583 HM583376 
Thylamys elegans MSB 87098 NK27606 HM583377 
Thylamys elegans MSB 133104 NK95111 HM583378 
Thylamys elegans MSB 209719 NK95354 KF164527 
Thylamys elegans MSB 209755 NK95436 KF164526 
Thylamys elegans MSB 209206 NK95622 KP994527 
Thylamys elegans MSB 209225 NK95677 KF164529 
Thylamys elegans MSB 209239 NK95691 KF164528 
Thylamys elegans MSB 133100 NK95971 KF164530 
Thylamys elegans MSB 210062 NK96571 KP994526 
Thylamys elegans MSB 133095 NK96763 HM583379 
Thylamys elegans MSB 133097 NK96791 HM583380 
Thylamys elegans MSB 210242 NK96879 KF164531 
Thylamys elegans SSUCMa519 KF164532 
Thylamys elegans SSUCMa520 KF164533 
Thylamys elegans UCK20 KP994524 
Thylamys elegans UCK21 KP994521 
Thylamys elegans UCK22 KP994531 
Thylamys elegans UCK23 KP994525 
Thylamys elegans UCK807 KP994529 
Thylamys elegans UCK809 KP994523 
Thylamys karimii MZUSP 32094 APC1561 HM583381 
Thylamys macrurus MZUSP 32094 APC932 HM583382 
Thylamys macrurus MSB 70700 NK27536 HM583383 
Thylamys pallidior MLP 24.X.01.3 AC47 HM583392 








APPENDIX 1 Continued 
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Species Museum Code Specimen Identifier Genbank Accession Number 
Thylamys pallidior DUS29 KF164522 
Thylamys pallidior EP434 KF164517 
Thylamys pallidior EP476 KF164515 
Thylamys pallidior HZP2691 KF164519 
Thylamys pallidior HZP3069 KF164511 
Thylamys pallidior HZP3079 KF164510 
Thylamys pallidior HZP3083 KF164509 
Thylamys pallidior JPJ1292 KF164524 
Thylamys pallidior MVZ115634 MVZ115634 HM583388 
Thylamys pallidior MVZ143696 MVZ143696 HM583389 
Thylamys pallidior MVZ173937 MVZ173937 HM583391 
Thylamys pallidior FMNH 162495 NBH7697 HM583384 
Thylamys pallidior MSB 57003 NK14721 HM583385 
Thylamys pallidior MSB 209881 NK96045 KF164513 
Thylamys pallidior MSB 210395 NK96067 KF164512 
Thylamys pallidior MSB 133108 NK96072 HM583387 
Thylamys pallidior OMNH 34908 OCGR7279 HM583406 
Thylamys pallidior OMNH 34911 OCGR7390 HM583408 
Thylamys pallidior PNG1055 KF164525 
Thylamys pallidior RRG34 KF164518 
Thylamys pallidior UCK481 KF164514 
Thylamys pallidior UP793 KF164521 
Thylamys pallidior VII485 KP994534 
Thylamys pusillus AMNH 260025- NK12539 HM583414 
Thylamys pusillus OMNH 23483 OCGRI1525 HM583369 
Thylamys sp. Giarla et al. 2010 MVZ116614 MVZ116614 HM583420 
Thylamys sp. Giarla et al. 2010 MVZ137896 MVZ137896 HM583423 
Thylamys sp. Giarla et al. 2010 HZP30 KF164554 
Thylamys sponsorius MSB 67015 NK23763 HM583441 
Thylamys sponsorius OMNH 34534 OCGR7432 HM583434 
Thylamys tatei MUSM23121 KF164555 
Thylamys tatei MUSM23253 KF164556 
Thylamys tatei MVZ135504 MVZ135504 HM583449 
Thylamys velutinus OMNH22284 OMNH22284 HM583450 
Thylamys velutinus OMNH37216 OMNH37216 HM583451 
Thylamys venustus MSB 63268 NK21664 HM583491 
Thylamys venustus AMNH 275428 NK22844 HM583468 


